Multi-resonant spinor dynamics in a Bose-Einstein condensate 
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We analyze the spinor dynamics of a 87 Rb F — 2 condensate initially prepared in the tuf = Zee- 
man sublevel. We show that this dynamics, characterized by the creation of correlated atomic pairs 
in mp — ±1, presents an intriguing multi-resonant magnetic field dependence induced by the trap 
inhomogeneity. This dependence is directly linked to the most unstable Bogoliubov spin excitations 
of the initial thf = condensate, showing that, in general, even a qualitative understanding of the 
pair creation efficiency in a spinor condensate requires a careful consideration of the confinement. 



Spinor Bose-Einstein condensates (BECs), consisting 
of atoms with non-zero spin, constitute an ideal scenario 
to investigate the interplay between internal and exter- 
nal degrees of freedom in a multi-component superfluid. 
The competition between spin-dependent collisional in- 
teractions, Zeeman effect, and inhomogeneous trapping 
results in an exciting range of fundamental phenomena. 
As a result, spinor BECs have attracted large attention 
since the pioneering experiments in optical traps [lj con- 
cerning their ground state properties @, 0, 0, la] and the 
spinor dynamics induced by the spin-changing collisions, 
which allow for a coherent transfer between different spin 
components [1,0]. 

Spinor BECs also provide exciting perspectives as 
novel sources of non-classical states of matter. In this 
sense, condensates initially prepared in the mp = Zee- 
man sublevel are especially fascinating 
that case, the creation of correlated pairs results in the 
growth of macroscopic populations in mp = ±1. This 
amplification process is ideally triggered by quantum spin 
fluctuations [lfj . Interestingly, it closely resembles para- 
metric amplification in optical parametric down conver- 
sion llj, opening exciting new routes for matter- wave 



squeezing and atomic Einstein-Podolsky-Rosen entangle- 
ment experiments 12, 13| . 

Correlated pair creation, and in general any spinor dy- 
namics, demands a significant rate of spin-changing col- 
lisions. In typical experiments these collisions are sup- 
pressed by the quadratic Zeeman effect (QZE) already 
in the presence of moderate magnetic fields [6[. How- 
ever, the influence of the QZE at low fields is far from 
trivial H Q H ES • In 

particular, spin-mixing can 
reach a pronounced maximum for low but finite fields. 
This resonance, contrary to those discussed below, has a 
non-linear character and has been explained in terms of 
phase matching flit ]. 

The understanding of the magnetic-field dependence of 
the pair creation efficiency is hence crucial for the char- 
acterization of novel spinor-based sources of non-classical 
matter waves. In this Letter we show that this depen- 
dence generally presents an intriguing non-monotonous 
character which is crucially determined by the trap in- 
homogeneity and cannot be explained from the physics 



of homogeneous BECs [18| . The pair creation efficiency 
is directly linked to the instability rate, which charac- 
terizes the exponential growth of the most unstable spin 
excitations of the initial BEC in mp = 0. This insta- 
bility rate presents pronounced maxima and minima as 
a function of the applied field, which result in the strik- 
ing multi-resonant magnetic field dependence of the pair 
creation efficiency (from mp — to mp = ±1) observed 
in our experiments. Along with these resonances we ob- 
serve characteristic magnetizati on p atterns which depend 
on both the magnetic field 0, [IJ [l9[ and the external 
confinement 1J)|, 2pj . 




FIG. 1: Illustration of the resonances in the spinor dynamics 
of a Bose-Einstein condensate initially in the \rriF — 0} state. 
For certain magnetic fields, the production efficiency of atoms 
in the \rriF = ±1} states is resonantly enhanced due to the 
excess energy released by the quadratic Zeeman effect. The 
column densities of all three Zeeman components in the plane 
of strongest confinement are shown in (a) and (c) in both a 3D 
and a density plot. Note that condensates in the |m_F = ±1} 
states are created in characteristic spatial modes depending 
on the particular resonance. 

To investigate this magnetic field dependence experi- 
mentally, we initially prepare a sample of 10 6 87 Rb-atoms 
in the F — 2,mp — 2 state in a crossed beam optical 
dipole trap at a wavelength of 1064 nm (see Refs. 2ll.l22| 
for a more detailed description). In the following, the 
mp = X states are denoted by \X). The spin polariza- 
tion of the atoms is maintained by a small guiding mag- 
netic field. The trap depth is lowered in multiple steps 
to evaporatively cool the atoms from a temperature of 
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FIG. 2: (a) Fraction of atoms transferred into \mp — ±1} 
within 21 ms as a function of the quadratic Zeeman energy q. 
To obtain this fraction, the number of transferred atoms per 
state was divided by the sum of transferred \m,F = ±1) and 
condensed \m,F = 0} atoms. Due to strong shot-to-shot fluc- 
tuations 15 independent realizations were averaged at each 
magnetic field. The error bars indicate statistical uncertain- 
ties. The blue line is a triple Gaussian fit to guide the eye. 
(b) Instability rate, given by the imaginary part of the most 
unstable spin Bogoliubov mode, corresponding to the pair 
creation efficiency into |m_p = ±1) (solid line). The maximal 
instability rate for an effective homogeneous case (dashed line) 
lacks any resonant features. 



2 /xK to 90 nK, resulting in a final trap depth of 150 nK 
and trapping frequencies of (176, 132,46) Hz. Typically 
a BEC of 7 x 10 4 atoms with a thermal fraction below 
40% is formed. Subsequently, the magnetic field is raised 
to 78 G and an adiabatic radio frequency passage is used 
to transfer the atoms to the |0) state within 5 ms. At 
this stage, the population of other spin states is below 
the detection limit of « 500 atoms. After the transfer, 
the field is lowered to 7.9 G within 15 ms. During this 
ramp, a strong gradient of « 50 G/cm is applied to ex- 
pel residual atoms in the mf / components from the 
trap. Subsequent to this purification, the field is quickly 
lowered to a specific value between 0.12 and 2 G within 

3 ms. This magnetic field strength is calibrated by mi- 
crowave hyperfine transitions and has a width of 4 mG 
due to current noise, current drift and residual gradients. 
Magnetic field gradients are below 0.2 mG/cm in all di- 
rections. The BEC is held at the chosen magnetic field for 
an adjustable time to allow for spin changing collisions. 
Finally, the dipole trap is switched off to let the atomic 
cloud expand freely. During time-of- flight, a strong mag- 
netic field gradient is applied in the vertical direction to 
separate the mj? components. The number of atoms in 
all five vtif components is measured by standard absorp- 



tion imaging along the weak trap axis. A bimodal fit 
to the measured density profile of the mp = compo- 
nent yields the condensed fraction and the temperature 
of the cloud. Typical density profiles after time-of-flight 
are shown in Fig. [TJ 

To evaluate the onset of the spinor dynamics, we re- 
strict our investigation to short spin evolution times and 
small populations in the |±1) states. Figure[5](a) shows 
the magnetic field dependence of the fraction of atoms 
transferred to the |±1) states after an evolution time of 
21 ms. The magnetic field strength B is represented by 
the quadratic Zeeman energy q oc B 2 which is the rele- 
vant energy scale. As expected 0], no atoms in the |±2) 
states were detected at this time scale. The pair-creation 
efficiency shows a striking mult i- resonant g-dependence. 
We show below that this dependence is directly linked 
to the interplay between Zeeman energy, spin-changing 
collisions and the inhomogeneous confinement of the en- 
semble. 

The onset of the spinor dynamics is best under- 
stood by considering the spin excitation modes of the 
initial condensate [18| . Depending on the quadratic 
Zeeman energy q and the atomic interactions, these 
modes may become dynamically unstable, triggering 
an exponential population growth in mp =/= states, 
hence initiating the spinor dynamics. An analysis 
based on spin excitation modes allows for a good un- 
derstanding of recent experiments in F = 1 87 Rb 
condensates, including spin-texture formation after a 
quench (20l . [23|, and - in combination with dipolar in- 
teraction - the instability of externally induced helices 
24 , [25|, . The initial state may be represented by a 
spinor wave function *f?o(r) = ("0-2, ip-i, ipo, ipi, ip2) T = 
(0, 0, no(r) 1 / 2 , 0, 0) T . The onset of the spinor dynam- 
ics (linear regime) is described by an operator H?(f,t) = 
(\&o(r) + <5\fr(r, i))e _vt , where // is the chemical poten- 
tial, 54?(r, t) — (Sif}-2,Siil}-i,SiilJo,5ipi,5ip2) T describes 
the spin fluctuations, and the density fulfils no(r) 3> 
(5ipj nF (f f )S^ mF (r)) . Up to second order in 8ip mF we ob- 
tain a Hamiltonian for the pair creation in |±1) [l^ |: 



J lm F =±r 

+ siefftf) [stlsjli + } . (1) 



Note that we may neglect the transfer into |±2) during 
the first stages of the dynamics due to the much lower 
transfer rate from |0). Hence the Hamiltonian is valid 
for the description of the initial pair creation for BECs 
in both F = 1 and F = 2. The effective Hamiltonian for 
the |±1) components is H e ff = —h 2 \/ 2 /2m + V e ff(r), 
with V e ff(r) = V(f) + (Uo + Ui)no(r) — //, where m is 
the atomic mass, and V(r) is the confining potential. The 
coupling coefficient £l e ff{f) = Uino(r) characterizes the 
pair creation induced by spin-changing collisions. The 
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interaction strengths are Uq = (7 go + IO52 + 18<74)/35, 
Ui = (-7g - 5g 2 + I2 54 )/35 for F = 2 and U = {go + 
252)73, U\ = (32 — <7o)/3 for F = 1 condensates. They 
are related to gp = A-KfPap /m, where is the s-wave 
scattering length for the collisional channel with total 
spin F (only even due to symmetry) 0. Note that for 
87 Rb Ui > for F = 2 whereas f7 a < for F = 1. 

For homogeneous BECs, V^// and fi e // are constants, 
and the spin excitations are plane waves with wavevector 
k and energy obtained after diagonalizing If. If 

the imaginary part Im(£fc(q)) is positive for some k, the 
BEC in |0) is dynamically unstable and pair production 
into the |±1) state occurs. By adapting the result of 
Ref. [18| to the general case of spinor BECs in F = 2 
and F = 1, we obtain three different regimes, which can 
be classified according to the value of q. BECs in the 
0} state are (i) stable for q > q cr + \q cr \ (where q cr — 
—Uino); (h) unstable for q cr < q < q cr + \q cr \, where the 
instability rate of the most unstable mode (with /c max = 
0) is A(q) = yq%. — (q — q C r) 2 /h; (hi) unstable for q < 
q cr with a constant instability rate A(g) = \q cr \/h for 
the most unstable modes h 2 k^ nax /2m — q cr — q. The q 
dependence of A(q) yields the magnetic field dependence 
of the pair creation efficiency. 

Although the homogeneous picture offers important in- 
sights, the observed g-dependence in trapped BECs is 
strikingly different, showing that the confinement must 
be considered to obtain even a qualitative understanding. 
This may be partially understood by comparison with the 
simple case of a F = 2 BEC in a box potential 1^(^ = 
for r < R, and V(f) = 00 otherwise. In the Thomas- 
Fermi regime the density no of the BEC in the |0) state is 
approximately constant for r < R, and the energies of the 
Bogoliubov modes are £^(g) — (e„ + q)(e n + q + 2[/ 1 n ), 
with e„ the energy of the n-th box state. Regimes (i) 
and (ii) are similar to the homogeneous case. However for 
regime (iii), occurring for 77(g) = \q\ — U\no > (eo + ei)/2, 
the most unstable mode is given by the level n whose en- 
ergy en is closest to 77(g). Contrary to the infinite system 
A(g) = \£n(q)\/h shows pronounced maxima and minima, 
reaching its maximum U±no only when 77(g) is resonant 
with a box state. 

For the actual experimental conditions not only the 
finite confinement but also the inhomogeneity of the po- 
tential and the related inhomogeneous Thomas-Fermi 
density profile no(r) contribute to the mode instability. 
This is best understood by introducing the eigenfunc- 
tions of the effective Hamiltonian H e f f(f) n (f) = e n 4> n (r). 
Note that the modes <fi n play a similar role to the box 
states in the simplified discussion above. After projecting 
Hm F {r) = J2n 4>n{r)a n , mF1 we may re- write the Hamil- 



tonian ([T|) in the form: 

H = (e„ + q) oJ,,ra r a n,mF + 

n,?7iF — ±1 

^ A n ^ (al l l a} l , _ 1 + a n / t _ia ni ij (2) 

where A n , n ' — J d 3 r^ e ff(f^(j> n (r)*(p n '(r). Note that 
the inhomogeneity of V(r) and no(r) plays a double 
role. First, V e ff(f) presents a non-trivial mexican-hat- 
like form (for U\ > 0). Second, contrary to the simpli- 
fied box potential, there is a significant coupling (A, un r) 
between different effective trap levels, induced by the in- 
homogeneity of Q e f f (r ) . 

After diagonalizing H we obtain the spin Bogoliu- 
bov modes £ v (q). Figure [2] (b) shows the magnetic 
field dependence of the maximal instability rate A(g) = 
max„|Im(£„(g))|/ft, for the experimental parameters of 
Fig. [2] (a). In the unstable regime with low |g| we may 
approximate the maximal instability rate by A(g) ~ 
V 9cr ~ (q ~ Qcr) 2 /h with an effective q cr (~ —30 Hz 
in figure [2|). However, as expected from the discussion 
above, this growth is not followed by a constant instabil- 
ity rate for larger \q\. On the contrary, A(g) shows pro- 
nounced maxima and minima, which lead to a strongly 
enhanced or reduced pair creation efficiency into |±1) 
due to the exponential nature of the growth. Note that 
the position of the maxima and minima of the observed 
pair creation efficiency in our experiments is in excellent 
agreement with the position of the calculated maxima 
and minima of the instability rate in Fig. [2] However, 
the evaluation of the resonance strength is more subtle. 
The experimentally observed resonance at low |g| seems 
to be stronger than the one at high \q\, whereas the calcu- 
lated instability rate is slightly smaller for the resonance 
at low |g|. This is not a disagreement, since the absolute 
population is not only determined by the growth rate but 
also by the initial condition. In fact, the resonant growth 
can be initiated by spuriously produced atoms in |±1) 
and by vacuum spin fluctuations. Since these triggering 
mechanisms contribute differently for the two resonances, 
they can fully explain the experimental results [Icil | . 

The agreement in the resonance position is maintained 
for different experimental parameters. In particular, we 
vary the atom number in our experiments by introducing 
an additional hold time of up to 30 ms before the purifi- 
cation step, which allows for a controlled loss of atoms. 
By varying this hold time, the total atom number N is 
changed without changing any of the other experimental 
parameters. Figure [3] shows the position of the resonance 
peak at low |g| (obtained from a Gaussian fit) for differ- 
ent numbers of atoms. The resonant g value increases 
with the number of atoms according to j /\r°- 36±0 - 02 . This 
exponent of w 2/5 closely resembles the dependence of 
the peak density on the number of atoms in the |0) state 
in the Thomas-Fermi regime. Hence the position of the 



4 




2.0 2.5 3.0 

Number of atoms / 10 4 

FIG. 3: Position of the lowest energy spinor dynamics res- 
onance, recorded for various atom numbers N of the initial 
condensate in the \rriF = 0) state. The blue line represents a 
fit to a power law TV 7 , yielding an exponent of 7 = 0.36±0.02. 



resonance at low \q\ depends linearly on the BEC density, 
as expected from the linear scaling of V e f / and il e f t with 
the density. 

Our results show that the pair-creation efficiency di- 
rectly reflects the instability of the spin Bogoliubov 
modes. This picture is general and may be applied to 
other experiments on spinor BECs. In particular, Fig. [4] 
shows the instability rate for the case of recent exper- 
iments with 87 Rb in F = 1 [9]. Our theory predicts 
a single resonance with a constant instability rate for 
< q < 6 Hz, and a clear decay of the pair-creation effi- 
ciency for q < (which was unexpected from the simpli- 
fied homogeneous approximation) , in excellent agreement 
with the reported results [9J]. We stress that the single- 
resonant character (as in Ref. @) or multi-resonant char- 
acter (as in our case) of the dynamics largely depends on 
atomic density and external confinement. 



trapping potential and the density of the BEC in the |0) 
state. We have shown that this inhomogeneity is cru- 
cial even for the qualitative understanding of the spinor 
dynamics. In particular, the trap-dependent modulation 
of the most unstable spin Bogoliubov mode is directly 
reflected by a strong enhancement or suppression of the 
spinor dynamics for particular magnetic field values. On 
the resonances, we observe the spontaneous emergence 
of characteristic magnetization patterns [28]. It is also 
important to stress that, in particular for Rb BECs in 
F = 1, the magnetic dipole-dipole interaction may play 
a significant role not only in the spatial magnetization 
pattern 0, H3 but also in the pair production effi- 
ciency. 

The characterization of the pair creation efficiency con- 
stitutes a comprehensive description of a parametric am- 
plifier for matter waves. The next challenge is the eval- 
uation of the amplifier input which can be of classical or 
quantum nature. A realization where the presented para- 
metric amplifier acts on pure vacuum spin fluctuations 
constitutes a promising source of non-classical states of 
matter [loj ]. This will open exciting perspectives for 
the detailed analysis of matter-wave entanglement and 
squee zing during the parametric amplification of the |±1) 
pairs EH , possibly allowing for the_ production of 
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atomic Einstein-Podolsky- Rosen pairs 
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FIG. 4: Instability rate as a function of the quadratic Zeeman 
energy for the parameters of Ref. [jj, (uj x , w y , u z ) = 2n x 
(440, 39, 4.2) Hz and N = 2 x 10 6 atoms. A constant maximal 
instability rate for < q < 6 Hz is followed by a significant 
decay of the conversion efficiency for q < 0. 

In conclusion, the pair creation efficiency in a spinor 
BEC initially in the |0) state is characterized by a multi- 
resonant magnetic field dependence, which results from 
the non-trivial interplay between QZE, spin-changing col- 
lisions and the finite size and inhomogeneity of both the 
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